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FffiLD OF THE INVENTION 

[0001] The present invention relates to methods for analyzing three-dimensional brain images 
described by brain volume data and, more particularly, to methods for determining a location for 
a mid-sagittal plane in a three-dimensional brain image. 

BACKGROUND 

[0002] To study an organ, such as a human brain, a three-dimensional image, or 'Volume 
scan," may be taken. An exemplary volume scan may be taken using Magnetic Resonance (MR) 
Imaging. The result of a volume scan may be considered to be volume data. For a review of tiie 
basics of MR Imaging, see Joseph P. Homak, The Basics of MRI, 1997 (available at 
www.cis.ritedu/htbooks/mri), the contents of v^ch are hereby incorporated herein by reference. 
In particular, see the section of Chapter 12 entitled Volume Imaging (3-D Imaging). For a 
standard database of MR volume data examples, see www.bic.mni.mcgill.ca/brainweb. The 
volume data may be considered to be organized in slices. Three types of slices are typically 
considered particularly useful, including axial slices taken normal to and along a vertical axis in 
an axial plane, which divides the brain into top and bottom portions. The slices also include 
coronal slices taken normal to and along a longitudinal axis in a coronal plane, which divides the 
brain into anterior and posterior portions. The slices further include sagittal slices taken normal to 
and along a transverse axis in a sagittal plane, which divides the brain into left and right portions. 
Typically, each slice is considered to be a two dimensional array of pixel intensity values. 
However, a quality other than intensity (e.g., hue, saturation, etc.) may be associated wifli a pixel. 
The analog of a pixel in three dimensions is called a voxel. A voxel may be considered to have a 
size in each of three directions. 

[0003] In human brain anatomy, two cerebral hemispheres may be identified as well as an 
interhemispheric fissure, which is a longitudinal furrow in the midline between the two cerebral 
hemispheres. A plane, called the '"mid-sagittal plane," may be defined as the sagittal plane 
passing through the interhemispheric fissure of the brain. The mid-sagittal plane typically 
contains less cerebral structure than a sagittal plane through any part of the rest of the brain in the 
vicinity of interhemispheric fissure. 

[0004] When a volume scan of a given brain is taken, the position of the given brain within 
the co-ordinate system of the volume data depends on the position of the patient's head during 
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the acquisition of the volume scan image. It is known that the volume data may be analyzed to 
determine a location for the mid-sagittal plane for the brain volume data. Determining the 
location for the mid-sagittal plane allows for the subsequent determination of a function to re- 
orient the mid-sagittal plane to a predetermined location. The re-orientation function, which may 
include translation and rotation, may then be applied to all of the brain volume data to orient the 
brain volume data in a preset co-ordinate system for further analysis. 

[0005] Previously developed methods for determining the location for the mid-sagittal plane 
are based on the extraction of symmetry lines in axial or coronal slices of the volume data and 
use a priori information about slice direction. The basic classes of methods include a class of 
methods that are based on the interhemispheric fissure and a class of methods that are based on a 
symmetry criterion. 

[0006] The basic hypotheses imderlying the class of methods that are based on the 
interhemispheric fissure include a hypothesis that.the interhemispheric fissure of the brain is 
roughly planar and a hypothesis that the interhemispheric fissure provides a good landmark for 
further volumetric symmetry analysis. In this class of methods generally, the fissure is identified 
as a segmented curve in MR images. Then, a three dimensional plane is found using an 
orthogonal regression from a set of control pomts representing the segmented curve. 

[0007] The theory behind the methods based on a symmetry criterion is that the mid-sagittal 
plane maximizes the similarity between the brain image and its reflection, i.e., that the mid- 
sagittal plane is the sagittal plane with respect to which the brain exhibits maximum symmetry. 
Most of the methods based on symmetry share a common general scheme. First, an adequate 
parameterization is chosen to characterize any plane of the three-dimensional Euclidian space by 
a vector composed of a few coefficients. For each selected plane in a set of possible planes, an 
adapted similarity meastire (synunetry criterion) is determined for the original three-dimensional 
brain image and a three-dimensional reflection image, where the reflection image is determined 
with respect to the selected plane. Then, the set of possible planes is searched to find the plane 
having the maximum adapted similarity measure. The chosen symmetry criterion is often the 
cross correlation be^veen the intensities of voxels in the two three-dimensional images. 

[0008] Unfortunately, known methods may include such drawbacks as a requirement for some 
a priori information, such as whether a given set of sUces is oriented in an axial plane, a coronal 
plane or a sagittal plane. Additionally, the methods may be considered time consuming, 
orientation dependent and limited to small tilts in the data. 
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[0009] Clearly there exists a need for a method of detennining the location of a mid-sagittal 
plane in a three-dimensional brain image that obviates a requirement for a priori information and 
overcomes other shortfalls of the previously developed methods. 

SUMMARY 

[00 1 0] Initially, volumes of interest are defined within a three-dimensional brain image 
described by brain volume data. Measures (entropy or energy) are then determined for the slices 
in the volumes of interest associated with each of three directions. One of the three directions is 
identified as the sagittal direction and a slice of &e volume associated with the sagittal direction 
and having the optimal measure may be used to define a first estimate of the mid-sagittal plane. 
The first estimate of the mid-sagittal plane may then be used in an optimization technique to find 
a final estimate of the mid-sagittal plane. 

[001 1] In accordance with an aspect of the present invention there is provided a method. The 
method includes detennining a candidate sagittal direction for a brain image, the brain image 
defined by brain volume data in a three-dimensional space associated with first, second and third 
directions, the first, second and third directions being orthogonal to each other, the candidate 
sagittal direction being the closest direction of the first, second and third directions to an actual 
sagittal direction. The determining the candidate sagittal direction includes defining a first three- 
dimensional volume of interest of the brain image, obtaining brain volume data in the first 
volume of interest for a first plurality of slices in the first direction, defining a second three- 
dimensional volume of interest of the brain image, obtaining brain volume data in the second 
volume of interest for a second plurality of slices in the second direction, defining a third three- 
dimensional volume of interest of the brain image, obtaining brain volume data in the third 
volume of interest for a third plurality of slices in the third direction, determining a measure for 
each slice of the first, second and third plurality of slices, for each of the first, second and third 
plurality of slices, plotting the measure for each slice of the plurality of slices, sequentially along 
a plot axis corresponding to respective first, second and third directions, to produce first, second 
and Mrd measure plots and detennining, from the first, second and third measure plots, which of 
the first, second and third directions is the candidate sagittal direction with a candidate plurality 
of slices associated therewith. 

[001 2] In accordance with another aspect of the present invention there is provided a method 
of determining an approximate location for a mid-sagittal slice firom a plurality of slices of bram 
volume data. The method includes obtaining the brain volume data for the plurality of slices, 
where the plurality of slices is generally oriented parallel to a sagittal plane and the plurality of 
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slices is positioned along an axis that is normal to the sagittal plane, detennimng a measure for 
each slice of the plurality of slices, wherein each measure is determined based on the brain 
volume data contained in each slice and selecting a candidate mid-sagittal slice among the 
plurality of slices, based on identifying an optimal measure amongst the measures determined for 
each slice of the plurality of slices. 

[001 3] In accordance with a further aspect of the present invention there is provided a method. 
The method includes determining a candidate sagittal direction for a brain image, the brain image 
defined by brain volume data in a three-dimensional space associated wiHi first, second and third 
directions, the first second and third directions being orthogonal to each other, the candidate 
sagittal direction being the closest direction of the first, second and tfakd directions to an actual 
sagittal direction. The determining of the candidate sagittal direction includes obtaining brain 
volume data for a first plurality of slices in the first direction, obtaining brain volume data for a 
second plurality of slices in the second direction, obtaining brain volume data in for a third 
plurality of slices m the third dkection, determining a measure for each slice of the first second 
and third plurality of slices, for each of the first, second and third plurality of slices, plotting the 
measure for each slice of the plurality of slices, sequentially along an axis corresponding to 
respective first, second and third directions to produce first, second and third measure plots and 
determining, fi"om the first, second and third measure plots, which of the first, second and third 
directions is the candidate sagittal direction with a candidate plurality of slices associated 
therewith. 

[0014] Other aspects and features of the present invention will become apparent to those of 
ordinary skill in the art upon review of the following description of specific embodiments of the 
invention in conjunction with the accompanying figures. 

BRIEF DESCRIPTION OF THE DRAWINGS 

[001 S] In the figures which illustrate example embodiments of this invention: 

[0016] FIG. 1 illustrates steps in a mid-sagittal plane-determining method according to an 
embodiment of the present invention; 

[0017] FIG. 2 illustrates three plots, each plot showing a normalized energy measure in each 
slice of a set of slices in a volume of interest associated with one of three directions; 

[001 8] FIG. 3 illustrates a normalized plot of measures for exemplary slices m a volume of 
interest associated with the sagittal direction along wifii images of seven of the slices; 
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[0019] FIG. 4A illustrates plots of various measures for sagittal slices of a brain volume scan; 



[0020] FIG. 4B illustrates plots of various measures for sagittal slices of a brain volume scan; 
[0021] FIG. 5A illustrates an initial mid-sagittal slice; 

[0022] FIG. 5B illustrates a first final mid-sagittal slice based on a final mid-sagittal plane 
obtained using an entropy (Kullback -Leibler) measure and Nelder-Mead optimization using an 
initial mid-sagittal plane defined based on the initial slice of FIG. 5A for building a starting 
simplex according to an embodiment of the present invention; 

[0023] FIG. SC illustrates a second final mid-sagittal slice based on a final mid-sagittal plane 
obtained using an energy measure and Nelder-Mead optimization using an initial mid-sagittal 
plane defined based on thetinitial slice of FIG. S A for building a startii^ simplex according to an 
embodiment of the present inventioi^ 

[0024] FIG. 6A illustrates an initial mid-sagittal slice; 

[0025] FIG. 6B illustrates a first final mid-sagittal slice based on a final mid-sagittal plane 
obtained using an entropy (Kullback -Leibler) measure and Nelder-Mead optimization using an 
initial mid-sagittal place defined based on the initial slice of FIG. 6 A for building a starting 
simplex according to an embodiment of the present invention; 

[0026] FIG. 6C illustrates a second final mid-sagittal slice based on a final mid-sagittal plane 
obtained using an energy measure and Nelder-Mead optimization using an initial mid-sagittal 
plane defined based on the initial slice of FIG. 6 A for building a starting simplex according to an 
embodiment of the present invention; 

[0027] FIG. 7A illxistrates a sagittal slice that is affected by statistical noise; 

[0028] FIG. 7B illustrates a sagittal slice that is affected by radio firequency iohomogeneity 
noise; 

[0029] FIG. 7C illustrates a sagittal slice that is affected by partial volume averaging noise; 

[0030] FIG. 8A illxistrates a linear representation of a mid-sagittal plane on an axial slice of 
brain volume data, where the mid-sagittal plane has been obtained using an energy measure 
according to an embodiment of the present invention; 
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[003 1] FIG. 8B illustrates a linear representation of a mid-sagittal plane on an axial slice of 
brain volume data, where the mid-sagittal plane has been obtained using an entropy (Kullback - 
Leibler) measure according to an embodiment of the present invention; 

[0032] FIG. 9A illustrates a linear representation of a mid-sagittal plane on a coronal slice of 
brain volume data, where the mid-sagittal plane has been obtained using an energy measure 
according to an embodiment of the present invention; and 

[0033] FIG. 9B illustrates a linear representation of a mid-sagittal plane on a coronal slice of 
brain volume data, where the mid-sagittal plane has been obtained using an entropy (Kullback - 
Leibler) measure according to an embodiment of the present invention. 

DETAILED DESCRIPTION 

[0034] In overview, volumes of interest are defined within a three-dimensional brain image by 
slices of brain volume data associated with three orthogonal directions. A measure, which may, 
for instance, be an energy or an entropy measure, is determined for each of the slices. The one 
direction, among the three directions associated with the volumes of interest, that most closely 
corresponds to the sagittal direction is then identified. The slice, among the slices in the volume 
of interest associated with the identified sagittal direction, having the optimal measure is used to 
define a first estunate of the mid-sagittal plane. The first estimate of the mid-sagittal plane may 
then be used as input to an optimization technique, which operates until a convergence criterion 
is satisfied, at which point a final estimate of the mid-sagittal plane may be produced. 

[0035] FIG. 1 illustrates steps in an exemplary mid-sagittal plane-determining method. 
Initially, brain volume data is read into a processor (step 102). The brain volume data may be 
MR data or any other image data. Notably, MR data is typically provided in slices oriented in 
three orthogonal planes or axes, which may, for example, be labeled X, Y and Z. However, for 
the purposes of the methods described herein, if the image data is not provided in slices, then the 
image data may be processed to rearrange the image data into slices. As should be clear to a 
person of ordinary skill in the art, three dimensional volume data contams a set of two- 
dimensional slices normal to and along each of three orthogonal directions (X, Y, Z). This is a 
basic principle of tomographic imaging, use of which allows one to obtain a two dimensional 
slice of an object described by three dimensional data. Each slice normal to a given axis may be 
associated with a coordinate in the given axis. For example, a slice normal to the X direction may 
be associated with x = 5, since, for each pixel in the slice, 
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[0036] Notably, the brain volume data can be either isotropic or anisotropic. That is, each 
voxel in the brain volume data may have a xmiform size along all three orthogonal axes or may 
exhibit different sizes along the three orthogonal axes. 

[0037] The number of slices normal to the X direction may be labeled Xmax, the number of 
slices normal to the Y direction may be labeled Ymax, and the number of slices normal to the Z 
direction may be labeled Zmax- Notably, there may be a different number of slices normal to each 
direction. Note that typically brain volume data is provided as an even number of slices in each 
direction. 

[003 8] The processor may then define a volume of interest associated with each direction 
(step 104). The volume of interest associated with a given dhrection is preferably centered about a 
mid slice of the slices normal to the given direction. An assumption is made in such a preference 
that the mid-s^ttal plane will pass through the vicinity of the centre of the im^e of the brain 
defined by the brain volume data and therefore will pass through all three volumes of interest. 
This assumption is expected to be true if the brain scan is of the whole of the brain. 

[0039] The three volumes of interest may be defined by first defining an index for the mid 
slices normal to the X, Y and Z directions as mid_x, mid_y and mid_z, respectively. The index 
defined for mid_x may, for instance, be jc = Xniax/2. The index defined for mid_y may, for 
instance, hQy = Yinax/2. The index defined for mid_z may, for instance, be z = Zn,ax/2. An index 
for an initial slice and a final slice normal to a given direction may then be determined based on 
the index for the mid slice normal to the given direction and a difference value: Init_x = mid_x - 
20 mm; Final_x = midjc + 20 mm; Initj/ = mid_y - 20 nun; Finaljy = mid_y + 20 mm; Init__z = 
mid_z - 20 mm; and Final_z = mid_z + 20 rom. The difference value may be considered to be 
converted into a number of slices, given the dimension of the voxels normal to the direction of 
interest. 

[0040] The three volumes of interest can then be defined using the initial slice and the final 
slice. The three volumes of interest associated with the three directions include: a volxmie of 
interest associated with the X direction (VOIx), whose candidate slices include only the slices 
normal to the X direction associated with the x-coordinates ranging from Init_x to Final_x; a 
volume of interest associated with the Y direction (VOIy), whose candidate slices include only 
the slices noraial to the Y direction associated with the jMioordinates ranging fi-om Init_y to 
Finally; and a volume of interest associated with the Z direction (VOI2), whose candidate slices 
include only the slices normal to the Z direction associated with the z-coordinates ranging fi:om 
Init_z to Final js. 
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[0041] As will be apparent to a person of ordinary skill in the art, the volumes of interest may 
be defined external to the processor and the processor may simply receive the volumes of interest 
as input. 

[0042] It should be noted tfiat a difference value of 40 mm (20 mm on either side of the mid 
slice) has been chosen empirically. The Applicants have taken into account the following. 
According to Duvemoy, H. M, The Human Brain: Surface, Blood Supply, and Three . 
Dimensional Sectional Anatomy, second edition, Springer, Vienna, 1999, p.302, the sagittal slice 
20 mm apart firom the medial sagittal slice has large manifestation of white matter, it cross-sects 
&e atriiun and the occipital horn of &e lateral ventricle and, additionally, has cerebrospinal fluid 
(CSF) in the meningeal space. The medial sagittal slice of brain may be shown (see Duvemoy, H. 
M, The Human Bram: Surface, Blood Supply, and Three Dimensional Sectional Anatomy, 
second edition, Springer, Vienna, 1999, p323) to contain a smaller amoimt of white matter 
structures and a larger manifestation of CSF (in the interhemispheric fissure). These differences 
in the manifestation of white matter and CSF around the interhemispheric fissure allow the 
location of the initial mid-sagittal slice. However, it should be recognized that the 40 mm 
difference value is provided as an example only and it should be clear that another suitable 
volume range may be chosen. 

[0043] For each volimae of interest, a measure (for example, entropy or energy) may be 
determined (step 106, FIG. 1) for each candidate slice. For a given candidate slice hoimal to the 
X direction, a first exemplary energy measure, /i, may be detennined fix)m: 

where xj^t is the intensity value of a pixel at the location (/, A) within the given candidate slice, 
where j represents the row index of the pixel and k represents the column index of the pixel. 

[0044] A second exemplary energy measure, /2, may be determined fix)m: 
[0045] A third exemplary energy measure, /3, may be detennined firom: 



[0046] A fourth exemplary energy measure, Zj, may be determined firom: 



f 
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[0047] A fifth exemplary energy measure, Is, may be determined firom: 

[0048] Linear combinations of the above energy measures may also be used. Although the 
above energy measures relate to candidate slices normal to the X direction, it should be clear that 
corresponding structure may be used for determining energy measures for candidate slices 
normal to the Y direction and the Z direction. These energy measures produce a single value for 
each candidate slice in the volume of mterest. 

[0049] Where / is an intensity value of a pixel in a given candidate slice, which pixel may 
have an intensity value selected from an intensity range from 0 to 255 inclusive, andp/ is the 
probability that a pixel randomly selected from among the pixels in the given candidate slice has* 
the intensity value /, a first exemplary entropy measure, 7i(p/), may be determined from: 

2SS 

which may be considered a measure of uncertainty of/?,. This measure is the negation of the well 
known "Shannon's entropj^' measure. 

[0050] Where pio is the probability that a pixel randomly selected from among the pixels in a 
reference slice has the intensity value /, a second exemplary entropy measure, IiiPnPio) , may 
be determined from: 

255 

/2(A. Ao) = Z^' 'Plo) 
0 

which may be considered a measure of the relative entropy of the distribution pt with respect to 
the distribution of p^, and is known as the 'Kullback-Leibler distance" between pi andp^o. 

[0051] The above entropy measures are determined by sumining over a range of intensity 
values. In the equations given, there are 256 intensity values. The reference sUce normal to each 
direction may be the corresponding initial slice. Like the energy measures, these entropy 
measures produce a single value for each slice. 
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[0052] Once measures, such as energy or entropy, have been determmed for each volume of 
mterest, the values of the measures may be normalized to take a value in the range from zero to 
one (step 108, FIG. 1). FIG. 2 illustrates a normalized plot of measures for slices in each volume 
of interest as having a horizontal plot axis and a vertical magnitude axis. In particular, FIG. 2 
illustrates a normalized plot of measures for slices normal to the X direction 202 (i.e., for VOIx, 
the X-axis is the plot axis), a normalized plot of measures for candidate slices normal to the Y 
direction 204 (i.e., for VOIy, Hhey-eods is the plot axis) and a normalized plot of measures for 
candidate shces in the Z direction 206 (i.e., for VOIz, the z-axis is the plot axis). The only 
normalized plot of measures to have a bell-shaped curve may be identified as having a plot axis 
corresponding to the sagittal direction. In particular, a processor may, for instance, ascertain the 
bell shape by recognizing that values of the measiire monotonically increase until a point or a set 
of points and then monotonically decrease. In the exemplary plots of FIG. 2, the normalized plot 
of measures for candidate slices normal to the X direction 202, has a bell-shaped curve and, as 
such, the X-axis may be identified as most closely corresponding to the sagittal direction. 
Notably, as defined herein, the sagittal direction is normal to each sagittal plane. 

[0053] Alternatively, a greatest degree of symmetry may be recognized in the normalized plot 
of measures for candidate shces normal to the X direction 202 when compared to the normalized 
plot of measures for candidate sUces normal to the Y direction 204 and the normalized plot of 
measures for candidate slices in the Z direction 206. That is, an axis orthogonal to the plot axis 
may be found, in the normalized plot of measures for candidate shces normal to the X direction 
202, about which the plot exhibits a greater degree of symmetry than either of the other two plots 
exhibit about any axis orthogonal to their respective plot axes. As such, the X direction may be 
identified as corresponding to the sagittal direction. 

[0054] The peak of the normalized plot of measures havmg a plot axis identified as 
corresponding to the sagittal dhection may be considered to be an optimal measure and to 
correspond to a candidate slice close to the mid-sagittal plane. In this way, a coarse estimate of 
the location of the mid-sagittal plane can be made from the three plots. 

[0055] In o&er embodiments the energy or entropy values are not normalized before being 
plotted. If a normalization function is used to normalize tiie plots, it should be clear that the 
normalization function should be an increasing function to preserve tiie maximum in the data so 
that an optimal measure can be determined. (An increasing function is one in which the bigger 
the value of the argument the bigger the value of that function.) 
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[0056] FIG. 3 illustrates another bell-shaped normalized plot of measures 302. (Notably, the 
normalized plot of measures 302 is distinct from the normalized plot of measures 202 of FIG. 2) 
The normalized plot of measures 302 is defined over 20 exemplary candidate slices normal to the 
sagittal direction along with images derived from seven of the 20 candidate slices. The middle 
seven of the 20 exemplary candidate slices are illustrated, as indicated by reference numeral 304, 
to illustrate changes in brain anatomy around the peak of the plot of measures 302. The candidate 
slice at which the peak of the bell-shaped normalized plot of measures 302 occurs may be 
selected as an initial mid-sagittal slice (step 110, FIG. 1). The initial mid-sagittal slice may be 
used to define a first estimate of a location for the mid-sagittal plane. The coordinate value of the 
initial mid-sagittal slice, in the sagittal direction, may be labeled i\r for convenience. 

[0057] FIG. 4A illustrates plots of three energy measures and two entropy measures for 37 
candidate slices of a first brain volume scan normal to a known sagittal direction. As can be seen 
from the plots, three of the measures produce a peak (or optimal measure) in the corresponding 
plot at candidate slice 21, with the remaining two measures producing peaks at candidate slices 
19 and 20 (i.e., iV= 21). 

[0058] FIG. 4B illustrates plots of ttiree energy measures and two entropy measures for 47 
candidate slices of a second bram volume scan (distinct firom the first brain volume scan of FIG. 
4A) normal to a known sagittal direction. As can be seen firom the plots, four of the measures 
produce a peak in the corresponding plot at candidate slice 27, with the remaining measure 
producing a peak at candidate slice 28 (i.e., iV= 27). 

[0059] Each optimal measure may, therefore, be considered to lead to the selection of an 
initial mid-sagittal slice, coordinates of pixels in which may provide a good startmg simplex for 
the optimization technique used to determine a location for a final mid-sagittal plane. 

[0060] In FIGS. 4A and 4B, a solid line represents Kullback-Leibler entropy measure, 
h(Pi>Pio) » ^ solid line connecting plus signs represents the fourth exemplary (log weighted) 
energy measure, /4, an alternating dot-dash line represents the third exemplary energy measure, 
/s, a dashed line represents Shamon entropy measure, /i(p/), and a dotted line represents a sixth 

r Y 

exemplary energy measure = ]^^yjk ~ S^y* • Applicants have achieved better results 



with the Kullback-Leibler entropy measure, l2(j>,»Pio) • However, it should be noted that the first 
exemplary energy measure, /i, is the least computationally intensive measure among the 
measures proposed herein. 
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[0061] It may be considered, for convemence, that the candidate slice selected as the initial 
mid-sagittal slice is normal to the X direction. The initial mid-sagittal slice may be considered to 
have four vertices, Pu Pi, Pz, Pa, which may be represented as: 

P,(x,y,2) = (Ml,l); 

P2ix,y,z)^(N,UZ„^; 

Pz(x,y,z)^(N,Ym^l);md 

PA{x,y,z)^{N,Y^7^. 

[0062] The initial mid-sagittal plane defined by the four vertices (Pi, P2, Ps^ P4) of the initial 
mid-sagittal slice may be used for building a starting simplex for an optimization process (step 
1 12, FIG. 1). The coordinates of the four vertices that define the initial mid-sagittal plane may be 
used as input for the optimization process. The optimization problem to be solved by the 
optimization process may be formed as an unconstrained optimization problem. Techniques that 
can be used to solve this type of optimization problem include the known grid-jfree algorithm of 
Nelder-Mead technique, the exhaustive search technique and the unconstrained optimization 
techniques available in Matlab™, a product marketed by The MathWorks of Natick, 
Massachusetts. 

[0063] As mentioned, an exemplary technique for solvmg an optimization problem that may 
be used to obtain a better estimate of a location for the mid-sagittal plane in a volume scan of a 
brain from the initial mid-sagittal plane is the Nelder-Mead optimization technique. The Nelder- 
Mead optimization technique is commonly used in non-linear regression problems and is knovm 
to operate under a number of rules, which may be expressed as steps. 

[0064] The Nelder-Mead optimization acts upon a simplex, which is a shape with m + 1 
vertices, where m is the number of parameters in the optimization problem. Thus, when the 
Nelder-Mead optimization is used to solve a two parameter optimization problem, the simplex 
has three vertices, i.e., the simplex is a triangle. The determining of a location for the mid-sagittal 
plane in a three-dimensional brain image may be considered to be a three parameter optimization 
problem (jn = 3). As such, the simplex for the Nelder-Mead optimization has four (m + 1) 
vertices. The initial mid-sagittal plane has four vertices. Pi, P2, P3, P4, as defined above, and may 
be used for building a starting simplex for a Nelder-Mead optimization. 
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[0065] Three of the four vertices of the initial mid-sagittal plane may be processed to 
determine a first vertex, 5b, in a three-dimensional parametric space. A plane in a three- 
dimensional parametric space is defined by three three-dimensional points. We may consider the 
initial mid-sagittal plane to be defined by vertices P2 and P3, Pr(xi, yi, Zf), Vertex Pi has a;; 
coordinate in common with vertex P2 and a z coordinate in common with vertex P3. As such, the 
first vertex, Sb> the three-dimensional parametric space may be defined by the set of three x 
coordinates for three vertices P/, P2, P3, that is, iSb = {xi, X2, Additionally, three more vertices 
in the three-dimensional parametric space may be defined by deviating firom iSb in each of the 
coordinates, that is: 

Si = {xi + A,;^C2,X3}; 

5i = {xi,X2 + A,X3};and 

[0066] From our exemplary vertices above, the four vertices in the three-dimensional 
parametric space are: 

52={^,iV+A,iV};and 
53 = {MMJ^+A}. 

[0067] To execute the Nelder-Mead optimization, a value for an objective function may be 
determined at each vertex of the simplex. The objective fimction may be, for instance, based on 
an original function such used to determine which slice of the volimie of interest should be 
selected as the initial mid-sagittal plane. However, while the Nelder-Mead optimization is 
designed to find a minimum, the original function used to select the initial mid-sagittal plane is a 
maximum for the slice closest to the mid-sagittal plane. As such, the objective function may be 
formed by negating the original fimction so that a point in the parametric space associated with 
the minimum value of the objective function is actually associated with the maximimi value of 
the original function. 
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[0068] For instance, where the original function is lu the objective function may be -/i, and 
determining the objective function for Si requires determining the sum of the intensities of the 
pixels at vertices: 

Pi(xi + A,;;,z) = (iV+A,l,l); 

i'2(x2, z) = (AT, 1, Zm^; and 

[0069] The Nelder-Mead optimization steps include: 

1) Reflect the vertex with the maximum value of the objective function through the centroid 
(center) of the simplex to locate a new vertex; 

2) If the value of the objective function at the new vertex is the minimum value among the 
values of the objective function at the other three vertices, expand the simplex and reflect 
further; 

3) If the value of the objective function at the new vertex is merely good, start at the top and 
reflect again; and 

4) If the value of the objective function at the new vertex is the maximum value among the 
values of the objective function at the other three vertices, compress the simplex and 
reflect closer. 

[0070] In the case of a three-dimensional parametric space, the simplex is a tetrahedron. The 
centroid may be found by averaging each of the coordinates of the three vertices that are not the 
vertex with the maximum value. 

[0071] The steps are repeated until a predetermined convergence criterion is satisfied, at 
which time a final mid-sagittal plane may be considered to have been defined. In particular, the 
best vertex in the three-dimensional parametric space defines three vertices in the three- 
dunensional brain volume scan space. Using these vertices, a plane equation (Ax-^By + Cz + D- 
0) may be formulated. Using this plane equation, an image of the optimal mid-sagittal plane may 
be constructed fiom the volmne data. The image of the optunal mid-sagittal plane may then be 
output. 



[0072] An exemplary convergence criterion can be defined as 
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a* ^ 

if=0 w 



where 

_ 2/« 

is the mean of the objective functions and^ is the value of the objective function at the 
vertex, S„. The Nelder-Mead optimization may be considered to have converged when a < s, 
where 8 is a predefined small number. 

[0073] The Nelder-Mead optimization may include rotation and translation to j5nd &e 
maximum value of the objective function. Notably, if the brain volume data is not rotated or 
translated relative to a pre-determined fiame of reference as provided, the translations and 
rotations associated with the optimizing step may not be required. However, such a case is 
expected to be rare. 

[0074] FIG. 5A illustrates an initial mid-sagittal slice 502. FIG. 5B illustrates a first final mid- 
sagittal slice 504 based on a final mid-sagittal plane obtained using the KuUback-Leibler measure 
and Nelder-Mead optimization using an initial mid-sagittal plane defined based on the initial 
slice 502 of FIG. 5A for building a starting simplex. FIG. 5C illustrates a second final mid- 
sagittal slice 506 based on a final mid-sagittal plane obtained using an energy measxire and 
Nelder-Mead optimization using an initial mid-sagittal plane defined based on the initial slice 
502 of FIG. 5A for building a startmg simplex. 

[0075] FIG. 6A illustrates an initial mid-sagittal slice 602. FIG. 6B illustrates a first final mid- 
sagittal slice 604 based on a final mid-sagittal plane obtained using the KuUback-Leibler measure 
and Nelder-Mead optimization using an initial mid-sagittal place defined based on the initial sUce 
602 of FIG. 6 A for building a starting simplex. FIG. 6C illustrates a second final mid-sagittal 
slice 606 based on a final mid-sagittal plane obtained using an energy measure and Nelder-Mead 
optimization using an initial mid-sagittal plane defined based on the initial slice 602 of FIG. 6A 
for buildmg a starting simplex. 

[0076] Once flie location for the final mid-sagittal plane has been determined, a function may 
be detennined to re-orient the final mid-sagittal plane to a predetermined location. The re- 
orientation function, which may include translation and rotation, may then be applied to all of the 
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brain volxime data using known techniques, such as the known AjBfine Transformation, to orient 
the brain volume data in a preset co-ordinate system for further analysis. 

[0077] The procedure of obtaining Magnetic Resonance Imaging (MBS) brain volume scan is 
known to be inherently marred by various types of noise. One type of noise is a statistical noise, 
which is random and of high frequency. Statistical noise increases as the signal strength increases 
or as the acquisition time is decreased. Another source of statistical noise is the electronics of the 
MRI system and can be attributed to thermal effects. A sagittal slice 702 that is affected by 
statistical noise is illustrated in FIG. 7A. Another type of noise, radio frequency (RF) 
inhomogeneity noise, may be attributed to the physical structure of the coils of the MRI system: 
RF inhomogeneity noise may be seen to occur due to many physical factors affecting tiie 
magnetization field generated by the MRI system. A sagittal slice 704 that is affected by RF 
inhomogeneity noise is illustrated in FIG. 7B. Partial volume averaging noise becomes more 
dominant as slice thickness increases. A sagittal slice 706 that is affected by partial voliraie 
averaging noise is illustrated in FIG. 7C. 

[0078] Advantageously, aspects of the methods of determining a location for a mid-sagittal 
plane in a three-dimensional brain image described above may be considered to function 
satisfactorily despite the presence of such noise as statistical noise, RF inhomogeneity noise and 
partial volume averaging noise. 

[0079] FIG. 8A illustrates a linear representation 802A of an optimized (final) mid-si^ttal 
plane on an axial slice of brain volume data 800, where the mid-sagittal plane has been obtained 
using the first exemplary energy measure, /i, outlined above. FIG. 8B illustrates a linear 
representation 802B of an optimized (final) mid-sagittal plane on the axial slice 800 of FIG. 8A, 
where the mid-sagittal plane has been obtained using the Kullback-Leibler entropy measure, 
hiPt»Pio) » oufii^ied above. 

[0080] FIG. 9 A illustrates a linear representation 902A of an optimized (final) mid-sagittal 
plane on a coronal slice of brain volume data 900, where the mid-sagittal plane has been obtained 
using the first exemplary energy measure, /i, outlined above. FIG. 9B illustrates a linear 
representation 902B of an optimized (final) mid-sagittal plane on the coronal slice 900 of FIG. 
9A, where the mid-sagittal plane has been obtained using the Kullback-Leibler entropy measure, 
IjiPiyPfo) 9 outlined above. 



[0081] Advantageously, aspects of the methods of determining a location for a mid-sagittal 
plane in a three-dimensional brain image described above obviate such drawbacks in previously 
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known methods, such as a requirement for a priori information on slice orientation. It may be 
considered that aspects of the present invention are orientation independent and are not limited to 
small tilts in the brain volume data. 

[0082] The foregoing describes the invention including preferred forms thereof. Alterations 
and modifications as will be obvious to those skilled in the art are intended to be incorporated in 
the scope hereof as defined by the accompanying claims. 



